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Summary. — 

Ten years have passed since the global solutions of advective accretion disks around 
black holes and neutron stars were first discovered. Since then they are enjoying 
support from observers almost on a daily basis, more so in recent days with the 
launching of very high resolution satellites. This review presents the development 
of the subject of advective accretion in last twenty five years leading to the global 
solutions and their applications. It also shows that apart from the standard Keple- 
rian disk features in most part of the accretion flow, future models must incorporate 
the essential features of the advective disks, such as the advection of energy and 
entropy by the flow, centrifugal barrier supported boundary layer of a black hole, 
steady and non-steady shocks, the bulk motion Comptonization of matter close to 
the black hole, outflows from the centrifugal barrier etc. Since black holes are 'black', 
methods of their identification must be indirect, and therefore, the solutions must 
be known very accurately. On the horizon, matter moves supersonically, but just 
before that it is subsonic due to centrifugal pressure dominated boundary layer or 
CENBOL where much of the infall energy is released and outflows are generated. 
In this review, we show that advective flow models treat accretion and winds onto 
black holes and neutron stars in the same footing. Similarly treated are the steady 
and time-dependent behaviour of the boundary layers of neutron stars and the black 
holes] Several observational results are presented which support the predictions of 
this advective accretion/outflow model. Signals from coalescing gravitating waves 
are affected by the sub-Keplerian flows as well. 

PACS 04.70.-S - Physics of black holes. 
PACS 97.10.Gz - Accretion and accretion disks. 
PACS 98.38.Fs - ets, outflows, and bipolar flows. 
PACS 33.20.Rm - X-ray spectra. 
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1. Introduction 

Standard accretion disk models of Shakura & Sunyaev [1] and Novikov & Thorne [2] 
have been very useful in interpretation of observations in binary systems and active galax- 
ies (e.g. [3, 4, 5]). The description of physical quantities in these models are expressed 
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analytically and they could be used directly. However, pressure and advection terms were 
not treated correctly, since the disk is terminated at the marginally sable orbit (three 
Schwarzschild radii for a non-rotating black hole). A second problem arose, when it was 
found out that the hot Keplerian disks of Shapiro, Lightman & Eardley [6] are viscously 
and thermally unstable close to the inner regions. Paczyhski and his collaborators [7, 8] 
attempted to include advection and pressure effects in the so-called transonic accretion 
disks, although no systematic study of global solutions was made. Global solutions of 
the so-called 'thick accretion disks' were possible to obtain only when the advection term 
is dropped (e.g., [9]). In these accretion disks, the flow is assumed to have practically 
constant angular momentum. Some exact solutions of fully general relativistic thick disks 
are discussed in Chakrabarti [10]. 

Meanwhile, observations of Cyg X-1 in the early seventies appears to indicate that it 
emits X-rays in two states [11, 12] which lead to speculations that probably the inner 
optically thin region of the disk emits the hard X-rays and the outer optically thick 
flows emit soft X-rays [13]. Ichimaru [14] correctly judging that the advection would 
be important close to a black hole, obtained a new optically thin solution (which in 
the self-similar limit arc known as Advection Dominated Accretion Flows of Narayan & 
Yi [15]) which includes heating, cooling and advection. This was done using a Newtonian 
potential. Ichimaru found that depending on the outer boundary condition, the solution 
could go over to the optically thick or the optically thin branch. There was no global 
solution, however. In the presence of a black hole, matter must be supersonic on the 
horizon, and, early attempts to find global solutions of viscous transonic flow (VTF) 
equations [8, 16] c;oncentrated much on the nature of the inner sonic point of these flows 
which is located around the marginally stable orbit. Meanwhile, in the case of inviscid 
adiabatic flow, an example of global solution was provided [17] where study of shocks 
similar to that in solar winds and galactic jets [18] was made. In the so-called 'slim-disk' 
model [19], using local solutions, it was shown that the instabilities near the inner edge 
could be removed by the addition of the advection term. This was done in optically 
thick limit. First satisfactory global solution of the governing equations in the optically 
thin or thick limit which include advection, viscosity, heating and cooling in the limit 
of isothermality condition was obtained by Chakrabarti [20, 21] where disk models of 
(single) temperature which become Keplerian far away were considered. 

In an earlier work [22] complete classification of global solutions of an inviscid, poly- 
tropic transonic flow (see, Fig. 4 of [22]) which showed that in some region of the 
parameter space, the flow will have multiple sonic points [23]. Within this region, there 
is a sub-class of solutions where Rankine-Hugoniot shock conditions are satisfied and 
standing shock waves are formed due to the centrifugal barrier. Four locations, namely, 
Xsi, (i ~ 1..4) were identified where these shocks could formally be possible, but it 
was pointed out that only Xs2 and x^s were important for accretion on black holes 
since the flow has to be supersonic on a black hole horizon and Xgi could also be im- 
portant for a neutron star accretion while .t_,4 was a purely formal shock location. In 
Chakrabarti [20, 21] viscosity was also added and complete global solutions in isother- 
mal VTFs with and without shocks, were found. In the language of Shakura-Sunyaev [1] 
viscosity parameter a, it was shown that if viscosity parameter is less than some critical 
value acr, the incoming flow may either have a continuous solution passing through the 
outer sonic point, or, it can have standing shock waves at Xs3 (following notations of 
Chakrabarti [22, 20, 21] if the flow allows such a solution in accretion. For a > acr, a 
standing shock wave at Xs2 persisted, but the flow now had two continuous solutions — 
one passed through the inner sonic point, and the other through the outer sonic point. 
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The one passing through the inner sonic point is clearly slowly moving in most of the 
regions, and therefore is optically thick. The one passing through the outer sonic point 
becomes optically thin at a large distance. A new branch of solutions which connected 
these two pieces is the standing shock wave solution which is partly optically thick (in the 
vicinity of the post-shock region) and partly optically thin (pre-shock flow coming out 
of a Keplerian disk). Later analytical and numerical works [24, 25, 26] showed that 
is stable, and that for a > acr the continuous solution passing through the inner sonic 
point is chosen. Most importantly, these solutions show that they could join (though not 
quite smoothly, since viscosity etc. were chosen to be constant) with a Keplerian disk at 
some distances, depending on viscosity and angular momentum [20]. 

Extensive numerical simulations of quasi-spherical, inviscid, adiabatic accretion flows [27, 
28, 30] show that shocks form very close to the location where vertically averaged model 
of adiabatic flows predict them [22] . The flow energy is conserved and the entropy gener- 
ated at the shock is totally advected into the black hole allowing the flow to pass through 
the inner sonic point. It was also found, exactly as predicted [22], that flows with positive 
energy and higher entropy form supersonic winds. In presence of viscosity also, very little 
energy radiates away (e.g.. Fig. 8 of [21]). Having satisfied with the stability of these 
solutions [24, 27, 25, 28, 29], a unified scheme of accretion disks was proposed [31, 25] 
which combines the physics of formation of sub-Keplerian disks with and without shock 
waves depending on viscosity parameters and angular momentum at the inner edge. The 
solutions remained equally valid for both the black hole and neutron star accretions as 
long as appropriate inner boundary conditions are employed. 

Subsequently, Chakrabarti [32] found that even when the isothermality condition is 
dropped, the flow topologies remained the same as in [20, 21]. Furthermore, it was 
shown that neutron star boundary layer could be studied as the post-shock region of the 
incoming flow itself (this was already pointed out [22]). The behaviour of the solution 
with viscosity was found to be intriguing. At very low viscosity or very high viscosities 
the flow does not have shocks, and sub-Keplerian, optically thin solutions emerge out 
of the Keplerian disks, very similar to the Ichimaru [14] solution. However, when the 
viscosity is intermediate, the solution passes through a standing shock for some region 
of the parameter space when 7 < 1.5 (a point completely ignored by workers who do 
not find these shocks thereby making the shocks 'controversial'). For other regions there 
are no shocks where abrupt changes in thermodynamical properties take place but more 
gradual changes are seen in the centrofugal barrier. 

Two dimensional solutions of advective disks were found in 1996 which showed the 
presence of funnels along the axis [33] . Such funnels were non-existent in the Advection 
dominated solution (see, Chakrabarti [34] for comparison of solutions.) Similar con- 
clusions, that the flow can't maintain its structure along the pole has been reached by 
Paczynski [35] and Park & Ostriker [36] as well. They showed that there is no way matter 
can be kept hot along the pole and matter must cool down. Exactly similar properties 
of the outflows from CENBOL has been used to explain the low frequency QPO of the 
black hole candiate GRS1915+105 [37, 38]. 

Several other developments were taking place in the theory of advective flows partic- 
ularly in the MHD limit and when the flow is non-axisymmetric. Chakrabarti [39] found 
all the global solutions of Weber-Devis type of equations. In this case, the nature of the 
magnetic field is predefined to be radial and azimuthal, but nevertheless the solutions 
indicated new types of sonic points in both accretion and winds. Instead of three sonic 
points (two X- types and one 'O' type), one obtains five sonic points (three X-type: slow 
magnetosonic, Alfvenic and fast magnetosonic and two 'O' type). Shocks were also found 
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analytically for the first time. Since then similar solutions have been found for cold flows 
where pressure and gravity effects were ignored (see, [40] and subsequent works of this 
group). However, till today, no new solution has been found other than those in C90c. 

When the advective flow is non-axisymmetric, spiral shocks may form in some regions 
of the parameter space [41, 42]. The later solutions are height integrated. Numerical 
simulations indicated that these self-similar shock solutions arc roughly correct [43]. 
Observationally, they may have been detected as well in M87. We shall discuss these 
observations later. 

2. — Milestones in Advective Accretion Flow Solutions 

Table I summarizes roughly the way the subject of advective accretion flows was 
developed in last twenty-five years or so. All self-similar solutions and those which are 
special cases of these solutions are excluded. In Table II, we show the summary of the 
observational evidence for the advective accretion flows. 

3. — Schematic pictures of accretion and outflow processes 

3'1. Building blocks of steady solutions. - With the advent of detailed and accurate ob- 
servations which demand more accurate solutions for proper explanations, it has become 
necessary to revise models of accretion and outflows. Figure 1 shows the classiflcation of 
all possible steady, inviscid solutions of transonic flows including schematic view of the 
disk around a black hole [22, 54]. This shows how an inviscid or weakly viscous flow look 
like. Beside each solution a schematic picture is given of accretion and outflows. Inward 
pointing arrows indicate the accretion solution and outward pointing arrows indicate the 
wind solutions. In (A) and (a), low energy and low angular momentum flow behaves 
like a spherical Bondi inflow and Parkar winds [75] well known in solar physics. In 
(B) and (b), a steady inflow or outflow solution does not have a shock but the inflow 
becomes hotter chie to slowing down at the centrifugal barrier. A time dependent simu- 
lation shows that an oscillating shock forms [30] which is accompanied by a non-steady 
outflow. In (C) and (c), the steady inflow solution has a standing shock and this has 
been verified by numerical simulations also [24, 27, 28]. Positive energy outflow is also 
found to be steady. In (D) and (d), the inflow passes through inner sonic point, but the 
outflow now has a steady shock [24]. In (E) and (e), neither the outflow nor the inflow 
has a shock in the steady solutions. But by analogy, it is expected that the outflow 
would have an oscillatory shock. This has not been tested yet. In (F) and (f), both the 
inflow and outflow have no shocks, and energetic steady solutions are present as both 
pass through the inner sonic point. In (G) and (H), incomplete solutions with closed 
topology are present. Thus no steady solution is possible. Simulations [76] indicated 
that flow is highly unsteady in these cases which are schematically shown in panels (g) 
and (h) respectively. Especially interesting is the solution of O*, whose disks could be 
'thick' not because of gas or radiation pressure, but because of turbulent pressure. 

In Figure 2, viscous transonic flow solutions are shown with constant viscosity pa- 
rameter [20, 21]. When viscosity is low, the shock, albeit weaker, still forms (see, panel 
A) and the disk joins with a Keplerian disk far out. The flow at the transition (Keple- 
rian to sub-Keplerian) radius need not be smooth (See, Fig. 3a for angular momentum 
distribution of the flow) and this may cause certain instabilities including oscillations 
and winds which carry away excess angular momentum [57, 77] . The nature of the disk 
corresponding to this is schematically shown in panel (a). In presence of higher viscosity, 
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Conclusions References 



Spectrum calculations using spherical flow [44, 45, 46] 

New optically thin branch of the solution [14] 

Multiplicity of sonic points [23] 

Systematic study of shocks in accretion [22, 20] 

Shocks in winds [20] 

The accretion shock is stable [24, 28] 

Classification of all inviscid accretion and wind solutions [22, 33] 

Positive Bernoulli constant required for shocks and winds [22] 

Global solution of viscous advective flows [20, 21, 32] 

Change of solution topology in presence of viscosity [21, 32] 

Fast- & slow magnetosonic points in MHD flows [39] 

Self-similar non-axisymmetric shocks [41, 42] 

Unification of Kepleriam and sub-Kepleriam flows under one scheme [31, 32] 

Post-shock flow has the similar properties as those of the hot ion tori [31] 

Gravitation waves emitted from binaries can be affected [47] 

Sub-Koplerian flows produce jets with more desirable properties [48] 

Vertically averaged flow model verified by two-dimensional simulations [27] 

Magnetic tension can cause disappearance of the inner disk [49, 50, 51] 

Necessity of a modified viscosity law in presence of advection [25] 

Keplerian/sub-Keplerian transition radius variation with viscosity [52, 32] 

Explanation of the spectral states by sub-Keplerian/Keplerian flows [52] 

Explanation of hard tail in soft state by Bulk motion Comptonization [52] 

Two domensional solution of advective disks including shocks [33] 

Self-consistent advective disk solutions in Kerr geometry [33, 54] 

Shocks may oscillate in a black hole accretion as in a White Dwarf accretion [53] 

Oscillating shocks can cause disappearance of the inner disk [30] 

Shock oscillation as the explanation of Quasi- Periodic Oscillations [53, 37] 

ADAF solution are special cases of advective disk model [32, 55, 56] 

Supersonic outflows are produced along funnel walls and at the transition radius [27, 57] 

Onflow production rate depends on compression ratio of the inflow [58, 34, 59] 

Neutron disks can be produced around a black hole [60] 

Nucleosynthesis can cause instabilities in advective disks [61] 

Bound matter of negative energy can be accelerated to infinity [62] 

Global Inflow-Outflow Solutions (GIOS) [63, 34, 64] 

Outflows can cause softening of spectra [65] 

Outflow rate depends on the spectral states [38] 

Outflows can cause quiescence states [64] 



the closed topology of Fig. l(B-C) opens up to accept matter straight from a hot flow 

provided the specific energy is positive. Better choice is to have the flow away from an 
equatorial plane where specific energy is likely to be positive. In case the flow comes 
from a cooler Keplerian disk, energy is negative, and the topology I* opens up to accept 
matter from it so that it may steadily enter the black hole through the inner sonic point 
(Fig. 2B). The corresponding disk is schematically shown in Fig. 2b. In Fig. 3b, the 
resulting angular momentum distribution is seen. If the viscosity is steadily reduced, the 
transition radius recedes outwards. It is to be noted that the viscosity parameter a we 
are discussing here is only that is valid in the sub-Keplerian flow. It has no relation to 
the a parameter in the Keplerian disk itself. 
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Table II. - Observational evidence for advective accretion flow solutions 



Conclusions 


References 


Hot radiation from accretion shock fits AGN spectra better 


[66] 


M87 Data [67] is best fitted with spiral shock solution 


[68] 


Timing analysis matches with two-component advective disks 


[69] 


Bulk motion Comptonization explains hard tails in soft states 


[70, 71, 72] 


Inner edge of the disk is around 10-20 Schwarzschild radii 


[82, 83, 84] 


Spectral evolution during novae outburst are explained 


[73] 


Quasi-Periodic Oscillation (QPO) is due to oscillation of shocks 


[37] 


Duration of QPO is related to frequency of QPO 


[37] 


Outflows lire produced witliiii iOOi?,, of tlie disk 


[74] 




Fig. 1: Nature of the solutions of inviscid advective disks and the corresponding 
flow behaviour. See text for details. 
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Fig. 2: Nature of the solutions of viscous advective disks and the corresponding 
flow behaviour: with shocks (A and a) and without shocks (B abd b). 



In realistic cases, neither the viscosity parameter is constant, nor the accretion rate 
remains constant. More importantly, one needs to consider processes such as Coulomb 
coupling between ions and electrons, bremsstrahlung, Comptonizaion etc. Also important 
is the availability of the driving forces (such as thermal pressure, centrifugal, magnetic 
etc.) to form outflows. Depending on these factors, a realistic flow would be made up of 
combinations of these basic building blocks of the inflow and outflow. 

One relavent matter in this regard is that ratio of the outflow rate to the inflow rate 
Rm depends on the strength of the shock, when it forms [58, 59]. One could argue that in 
soft states, when the accretion rate is high, the inner sub-Keplerian region is cooled down 
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Fig. 3: Nature of anguar momentum distribution (solid) in an advective disk 
close to a black hole as a function of the logarithmic radial distance. Two cases 
are same as (a) and (b) of Fig. 2. Dotted curves denote Keplerian distribution 
for comparison. 
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Fig. 4: Variation of the ratio of the outflow rate to the inflow rate Rm with the 
compression ratio 7? of the shock. It is believed that in soft states, 7? ~ f and 
thus the outflow rate is also very low. In hard states, shocks are stronger and 
compression ratio is higher. Here, Rm is also low. For intermediate J? ~ 3, Rm is 
highest, outer sonic point of the flow rc is closest to the black hole. In this case, 
the outflow may be periodically cooled and it may therefore be blobby [38, 37]. 
Dotted curve shows the actual outflow rate Mout when Mi„ varies as 1/R a 
plusible choice for an advective flow. 



[52] and the shock is as good as non-existent and there is no driving force to produce 
outflows [38]. In the hard states, on the contrary, when the accretion rate is low, the 
stronger shocks also causes very low outflow rates. In the absence of shocks, driving force 
is weak, and the outflow may be negligible (in the absence of strong magnetic flelds). 
For intermediate shock strength, the outflow rate is high and it could be cooled down 
by soft photons from the pre-shock Keplerian disk. Once the flow is cooled down, it 
would be super-sonic and fly away but the terminal velocity would be low. This process 
may take place repeatedly [38, 37] and cause low frequency variations of the light curve. 
In this so called 'flare' state, the outflow would be 'blobby'. The variation of Rm with 
shock compression ratio is shown (solid) in Fig. 4 with a suggestive dependence of Rm 
shown on the spectral state of the flow. The dashed curve shows the variation of the 
actual outflow rate when the accretion is assumed to vary as Mi„ = ^ . This was done in 
keeping with the fact that in soft states {R ^ 1) the inflow rate is around Aii„ ~ MEdd 
and in very hard states (i? ~ 7), Min ^ O.OlMEdd- 

3'2. Building blocks of non-steady solutions. - Two major types of oscillations of the 
centrifugal pressure supported boundary layer is discussed in the literature, and it is 
possible that both are important. Figure 5 shows the flow conflgurations at different 
times (measured in units of 2GM/c). This happens when the steady shock solution is 
absent even when the inner sonic point exists [30]. In this case, the flow searches for a 
steady solution by flrst generating entropy through turbulence and then trying to pass 
through the inner sonic point non-existent in a steady flow. Figure 6(a-b) shows the 
nature of the variation of the entire disk when bremsstrahlung cooling is added [53]. 
This is a prototypical case, when the steady shock solution exists but the cooling time 




Fig. 5: Oscillation of shocks and corresponding periodic changes in velocity 
field when flow parameters are 'wrong'. Sec [30] for details. 



in the post-shock/corona region roughly agrees with the infall time in this region. An 
oscillation of similar kind is seen at the transition radius between Keplerian and sub- 
Keplerian flows as some matter is removed as winds from this region [57, 77]. This is 
present even when no shocks are produced. A third type of oscillation may be set in 
due to the quasi-periodic coohng of the outflow (provided it is high as is the case for 
average shock strengths R ^ 2.5 — 3 see, Fig. 4) due to enhanced interception of the 
soft-photons emitted from a Keplerian disk. It can be assumed that the quasi-periodic 
oscillation (QPO) in X-rays observed in black hole candidates is the result of one or more 
of these different types of oscillations since other types of oscillations, such as those due 
to trapped oscillations [78] or diskoseismology [79] are incapable of modulating X-rays 
with large scale amplitude. 

3'3. Building blocks of magnetized aceretion and winds. - All possible solution topolo- 
gies of transonic flows in presence of magnetic field are in Chakrabarti [20, 39] and are 
not discussed here. The topologies would not change, but the actual solution would be 
modified slightly when Grad-Shafranov equations are solved. 
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Fig. 6a: Oscillation of shocks and corresponding periodic changes in velocity 

field when when flow parameters are such that the infall time matches with the 
cooling time. The corona is prominent at this phase. See [53] for details. 
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Fig. 6b: Same as in Fig. 6a but after T/2 time when the corona is collapsed. 

See [53] for details. 
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3'4. Combinations of building blocks. - One could combine results of Figs. 1-6 with 
to obtain realistic disk behaviour around a compact object. Up until 5-6 years ago, it 
was believed that the hard radiation is the result of Comptonization of soft photons 
by 'Compton Clouds' floating around the disk or by hot corona (Fig. 7). Chakrabarti 
and Titarchuk [52], for the first time pointed out that the so-called 'Compton Cloud' 
is the blottcd-out, sub-Keplerian inner edge of the Originally disk itself! (This is to be 
contrasted with the earlier models [14, 6] which considered only the Keplerian disks.) 
Thus, a new paradigm of accretion disks, based on actual solution of advective flows, 
started. Today, this picture is universally adopted in most of the models of accretion 
flows. One can have one of the several types of flows shown in Fig. 8. The so-called 
'hard' state is shown in Fig. 8(a-b) and the so-called 'soft' state is shown in Fig. Sc. 
The difference is that in the hard state, Keplerian flow rate coiild be very small while 
the sub-Keplerian rate could be very high [52]. In the soft state, it is the opposite. 
The de-segregation of matter into these two types of rates are believed to be due to the 
fact that the flow closer to the equatorial plane is likely to be more viscous and with 
larger Shakura-Sunyacv parameter a, and therefore is likely to be Keplerian. However, 
flows away from the equatorial plane may have smaller a and therefore they deviate 
from a Keplerian disk farther out [32]. Contribution to this sub-Kcplcrian flow may also 
have come from wind accretion. This sub-Keplerian flow may form shocks at around 
10 — 2{)Rg. In the soft state, shocks may be nominally present, but would be cooler due 
to Comptonization, and would be as good as non-existent. The combined picture in these 
two states was shown in Chakrabarti & Titarchuk [52]. As far as the wind and outflow 
formation goes, the result strongly depends on the compression ratio (i?) of the flow at 
the shock as discussed in Chakrabarti [34, 59] where an analytical expression (Fig. 4) for 
the ratio Rm of the outflow (assumed isothermal) to inflow rates is provided. When the 
compression ratio is unity, i.e., in the shock-free case (i.e., soft state), the outflow rate 
is negligible. Thus, in soft states, no outflow is expected. In the hard state, the shock 
is strong, and the outflow may be significant. In the intermediate state outflow rate is 
highest. In Das & Chakrabarti [64] it was shown that should generally depend on 
infiow rates as well. If a shock does not form and the accretion rate is low, the spectrum 
will be harder and no quasi-periodic oscillations (QPO) would be seen [73] in hard X- 
rays. This type of disks is shown schematically in Fig. 8a. If the shocks form, QPOs 
would be seen in hard X-rays, but softer X-rays, arising out of a Keplerian disk will 
not be participating in the QPOs. This type of disk is schematically shown in Fig. 8b. 
When the viscosity and fiow accretion rate is large, the sub-Keplerian region shrinks as 
it is cooled down by thermal Comptonization. This is schematically shown in Fig. 8c. 
Only the bulk motion Comptonization will take place in flows at around 1 — 3i?g thereby 
producing hard tails in soft states [52]. Outflow also softens the spectra mimicking those 
from bulk motion Comptonization [65], but as discussed above, it is unlikely that there 
would be strong outflows in the soft state due to the lack of driving forces. 

3'5. Blobby Jets?. - We now discuss the conditions under which the jets could be 
blobby. Such jets were predicted purely from theoretical point of view, in the context 
of anti-correlation found in the X-ray and radio fluxes in the black hole candidate GRS 
1915-1-105 (see, §4.7-4.8 of [55]). Another easy way is to use magnetic tension. When 
magnetic fleld is dragged in by the accretion disks, it could be sheared and amplifled in 
the post-shock region till an equipartition is reached. A strong field in a hot gas feels 
magnetic tension ('rubber-band effect') and is contracted catastrophically evacuating the 
post-shock flow, i.e., the inner part of the accretion disk. Quoting [55]: " Formation of 
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Fig. 7 : Cartoon diagram of a Keplerian disk with clouds of electrons floating on 
both sides. This was the conception of black hole accretion flow before advective 
disk solutions were applied. 



coronae in an accretion disk is not automatic, i.e., presence of magnetic field inside a disk 

docs not automatically imply that a hot magnetic coronae would form. The formation 
requires the ability of the accretion disks to anchor flux tubes inside the disk. This means 
that the disk should have an internal structure akin to the solar interior and the entropy 
must increase outwards. If the entropy condition is proper, the coronae would form, 
otherwise it would come out of the disk as a whole without causing any random flare. In 
former case, there would be sporadic flaring events on the disk surface as in the case of the 
sun, whereas in the latter case, the collapse of fields in the funnel would cause destruction 
of the inner part of the disk and formation of blobby radio jets. Detailed observation of 
GRS1915+105 shows these features [80]. Since the inner part of the accretion disk could 
literally disappear by this magnetic process, radio flares should accompany reduction of 
X-ray flux in this objects. Since the physical process is generic, such processes could also 
be responsible for the formation of jets in active galaxies and similar anti-correlation may 
be expected, though time delay effects are to be incorporated for a detailed modeling." 
In the absence of the magnetic field, periodic evacuation by oscillating shocks can also 
take place if the flow parameters are 'wrong' (see, [30] for details). This phenomenon 
may have been observed very convincingly in GRSf 915+105. In addition, the outflow 
may be periodically cooled. In Fig. 8d, we show disks in a possible flare-state as observed 
in GRS 1915+105. Here shock strength is intermediate and the outflow is blobby. 

4. Observations of Advective Accretion and Winds 

Several observations are now available which could be explained only with advective 
disks since normal Keplerian disks appear to be insufficient. 

4T. Sub-Keplerian motion on a large scale. - Since shocks are transitions from Su- 
personic to sub-sonic motion, and since supersonic flows are sub-Keplerian ( [32, 55]), 
any presence of shocks in a disk would indicate sub-Keplerian motions. Sub-Keplerian 
flows rotate slower, and velocity predicted from Doppler shifted disk emission lines would 
correspond to a higher Central mass. Chakrabarti [68] pointed out that the disk around 
M87 contained shocks and therefore the flow must be sub-Keplerian. As a result, the 
mass of the central object was found to be around 4 x 1O^M0 rather than 2 x IO^Mq 
as predicted by Harms et al. (see, [68] for references) with calculations purely based on 
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(a) 




HARD STATE-A 



Fig. 8a: Cartoon diagram of an advective disk when the building blocks of Fig. 
l(B-C) and Fig. 2(A-B) were used. A hard state is expected with or without 
QPOs. Wind may be formed at a low rate (Fig. 4). 




HARD STATE- B 



Fig. 8b: Another possibility of hard states where shocks are not formed in 
accretion. Building blocks of Fig. 1(A, D) and Fig. 2B were used. A hard state 
is expected without QPOs. Winds may be formed at a low rate (Fig. 4). 




Fig. 8c: Cartoon diagram of the advective disk with a soft state where shocks, 
if formed, are cooled down in accretion. Building blocks may be any of Fig. 
l(A-F) and Fig. 2B were used. A soft state is expected without any QPO. N 
significant wind is possible (Fig. 4). 



14 



S. K. CHAKRABARTI 



Soft x-rays Hard X-rays 




FLARE STATE 



Fig. 8d: Cartoon diagram of the advective disk in a flare state where shocks 

of intermediate strength may form but since the outflow rate is high (Fig. 4) , it 
may be periodically cooled to produce stange behaviour as in GRS1915-I-105. 



Kcplcrian motion. The fact that shock-ionization causes the emission processes on M87 
disk has been stressed recently by several others [81]. 

4'2. Sub-Keplerian motion on a small scale. - It is usual practice to assume that the 
inner edge of a Keplerian disk extends till 3Rg, the marginally stable orbit. However, 
the advective disk models [55, 32] show that the inner edge would be typically around 
~ 10 — 2QRg where the CENBOL should form. There are overwhelming evidence today 
that this is indeed the case [82, 84, 83]. This sharply contrasts predictions of ADAF 
model which often models the inner edge a distance of several tens of thousands (even 
at hundreds of thousands) of Schwarzschild radii! 

4'3. Power-law hard radiation in very high states. - Chakrabarti & Titarchuk [52] 
pointed out that when the accretion rate is relatively high, the electrons in the sub- 
Kcplcrian region become cooler and this region becomes practically indistinguishable 
from that of Keplerian disk. However, very close to the black hole horizon, matter moves 
with almost velocity of light and deposits its bulk momentum onto the photons thereby 
energizing these photons to very high energy forming a power-law. This power-law is the 
hall-mark of all the known black holes [70]. The success of this model crucially hangs 
on the transonic flow solution which utilizes the fact that the inner boundary condition 
is independent of the history of incoming matter. Another possiblity of producing hard 
spectra in soft states would be to use softening property of outflowing winds [65] but 
since outflow should be weaker in soft states (Fig. 4), this procedure is not obvious. 

4"4. Quasi-Periodic Oscillations from black hole candidates. - X-rays from galactic 
black hole candidates often show persistent oscillations which are quasi-periodic in nature. 
Advective disk solutions do allow oscillations (Figs. 5-6), especially in the X-ray emitting 
regions. The observation of these oscillations during the transition of states is the triumph 
of the advective disk model [85]. ADAF models, unlike advective disk solutions described 
above do not allow shock formation [86] and therefore cannot address the issues of quasi- 
periodic oscillations at all. 

4'5. Outflows and their effects on the spectral properties. - As discussed in §3.1 and 
§3.2, there are outgoing transonic solutions which are counterparts of transonic accretion 
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flows onto black holes [22]. These outflows may or may not have shocks, but they 
invariably have positive binding energy at a large distance. ADAF solution recently 
came to the same conclusion (that Bernoulli parameter should be positive [87]). However, 
when the outflow is accelerated by external agency such a conclusion is wrong as has been 
recently shown by Chattopadhyay and Chakrabarti [62] . They showed that even matter 
with negative initial binding energy could be pushed to infinity (large distance) when 
radiation momentum is systematically deposited on the flow! Recently, to overcome 
some of the shortcomings, ADAF has also tried to produce some self-similar outflows 
(ADIOS [88]) Being self-similar, jets are to come out all over the disk, as opposed to 
regions close to the black holes as predicted by our advective disk solutions. Recent high 
resolution observations of jets in M87 strongly suggest that they are produced within a 
few tens of Schwarzschild radii of the horizon [74]. strongly rejecting ADAF and ADIOS 
models for the outflows, which has no special length scale at these distances. Given 
an inflow rate, one is now capable of computing the outflow rate when the compression 
ratio at the shock surface is provided [63, 59]. This solution naturally predicts that the 
outflow must form at the CENBOL. The Globally complete Inflow-Outflow Solutions 
(GIOS) were also found [63]. In presence of winds, the spectra is modified: outflows take 
away matter and reduces the number of electrons from the post-shock (CENBOL) region 
and thus they are cooled more easily by the same number of soft photons emitted by the 
pre-shock Keplerian flow. This makes the spectra much softer [58] in presence of winds. 

4'6. Quiescence states of black holes. - Chakrabarti [34] and Das & Chakrabarti [64], 
pointed out that in some regions of the parameter space, the outflow could be so high 
that it evacuates the disk and forms what is known as the quiescence states of black 
holes. A well known example is the starving black hole at Sgr A* at our galactic centre 
whose mass is ~ 2.6 x IO^Mq and the accretion rate is supposed to be around ~ 10~^Mq 
yr~^ which is much smaller compared to the Eddington rate. Quiescence states are also 
seen in stellar mass black hole candidates such as A0620-00 and V404 Cygni. Another 
way of producing these states is to use well known viscous instability in an accretion disk 
as used in models of dwarf- novae outbursts [52]. 

4'7. On and Off-states during QPOs in black holes. - The black hole candidates 
GRS1915-I-105 displays a variety of behaviour: usual high frequency QPOs are frequently 
interrupted by low frequency oscillations. While the QPO frequency can be explained by 
the shock oscillations, the switching of on and off states is explained by the duration in 
which extended corona becomes optically thick. It is possible that outflows are slowed 
down by this process and the matter falls back to CENBOL, extending the duration of 
the 'on' state (which, if exists, is found to be comparable to the duration of the 'off' 
state [89, 90]). The QPO frequency does evolve during this time scale and it is suggested 
that this is due to the steady movement of the inner edge of the Keplerian accretion 
disk [91] or the steady movement of the shock itself in viscous time scale [92]. Generally 
speaking, shocks may show large scale oscillations [53, 30] but the wind plays a vital 
role in intercepting larger number of soft photons, characteristics of soft state. The 
correlation between the duration of the off-state and the QPO frequency has been found 
to agree with the observations [92]. 

4'8. Relationship of outflows and the black hole states. - Fig. 4 shows the possible 
relation between the outflow rates and the spectral states of the black holes. Recently, 
from observational point of view this picture has gained some support as well [93]. 
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4'9. Effect of advective disks on gravitational wave. - It was pointed out [47, 94] that 
since advective disks close to a black hole need not be Keplerian, it would affect the 
gravitational wave properties of a coalescing binary. The angular momentum loss to the 
disk by the companion in comparison to the angular momentum loss to gravitation wave 
was found to be very significant [94] . The variation of the signal is shown in [58] One of 
the exciting predictions of this scenario is that since the spectrum of an accretion disk 
contains a large number of informations (e.g. mass of the central black holes, distance 
of the black hole, accretion rate, and viscosity parameter) a simultaneous observation 
of the electromagnetic spectrum from the disk and the gravity wave spectrum (which 
also must depend on those parameters, except possibly the distance) should tighten the 
parameters very strongly. This issue is being investigated and will be reported shortly. 

Recently, ADAF model of extremely low accretion rate was used to repeat these 
computations, and not-surprisingly, no significant change in gravitational signal was 
found [86]. If little matter is accreting, it is as good as having no accretion disk at 
all. ADAF model ideally valid for zero accretion rate systems and not found to be true 
in any of the realistic accreting systems. As discussed earlier, evidences of the oscillating 
shocks, bulk motions etc. (which are absent in ADAF) are abundant. Thus it is likely 
that the gravitational signal would be affected exactly in a way computed by [94] 

5. — Concluding remarks 

When globally complete advective disk solutions were first introduced ten years ago 
which included centrifugal barrier and standing shocks, it was considered to be of aca- 
demic interest. Today, with the advent of accurate observational techniques, including 
those with high time resolutions, there is enormous evidence that these are more real- 
istic solutions of the flow. There is hardly a single paper in the literature today which 
does not use the advective disk solutions (sometimes not recognizing these point at all 
as in ADAF models.) As we have already noted, most of these improved observations 
cannot be explained at all using standard Keplerian disk or any simple variation of it. 
We believe that all future models of the disks must include the basic properties of this 
model including centrifugal barrier, standing shocks, outflowing winds, sub-Keplerian 
inner region etc. 

* * * 

The author thanks DST for a partial support for the project 'Analytical and numerical 
studies of astrophysical flows around black holes and neutron stars' and the organizers 
to partly supporting his trip to Rome and Pescara and covering the local hospitalities. 
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